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ABSTRACT 

I show that Tremaine- Weinberg (TW) measurements of bar pattern speeds are 
sensitive to errors in the position angle of the disc, PAdisc- I use an TV-body experiment 
to measure these errors; for typical random PAdisc errors, the resulting scatter in the 
measured values of the dimensionless bar speed parameter TZ (defined as the ratio of 
the corotation radius to the bar semi-major axis) is of order the observational. 

I also consider how the systematic PAdisc errors produced by disc ellipticities affect 
TW measurements. The scatter produced by these errors may be significant, depending 
on the ellipticity distribution. Conversely, by using the sample of TW observations, I 
find that an upper limit of the typical disc (density) ellipticity is 0.07 at the 90 per 
cent confidence level, which is in good agreement with previous measurements. 

Taken together, the random and systematic scatter suggest that the intrinsic dis- 
tribution of TZ of gas-poor early-type barred galaxies may be as narrow as that of the 
gas-rich later-types. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: kinematics and dy- 
namics — galaxies: structure — methods: observational 



1 INTRODUCTION 

Barred (SB) galaixies account for more than half of all high 
surface brightness disc galaxies (Knapen 1999; Eskridge et 
al. 2000). Recent observational and theoretical studies have 
focused on the pattern speed of bars, fip. The quantity of 
greatest interest is TZ = Dl/cib, where Dl is the corotation 
radius and as is the semi-major axis of the bar. A self- 
consistent bar must have TZ > 1 (Contopoulos 1980); bars 
with 1.0 < TZ < 1.4 are termed fast, while slow bars have 
larger TZ. Because bars have strong quadrupole moments, 
they lose angular momentum efficiently in the presence of 
a dense dark matter halo (Weinberg 1985), slowing down 
in the process; fast bars therefore have been interpreted as 
evidence for maximum discs (Debattista & Sellwood 1998, 
2000, but see also Valenzuela & Klypin 2002). Thus the ac- 
curate measurement of TZ in SB galaxies is of interest. 

Bar pattern speeds can be most reliably measured when 
kinematic data are available. One method relies on the de- 
pendence of the gas flow pattern on fip, particularly at the 
shocks in the bar region. Hydrodynamical simulations can 
therefore recover Qp; these find fast bars {e.g. van Albada 
& Sanders 1982; Athanassoula 1992; Lindblad & Kristen 
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1996; Lindblad et al. 1996; Weiner et al. 2001). An alter- 
native method, which measures flp directly, was developed 
by Tremaine & Weinberg (1984). Until now, the Tremaine- 
Weinberg (hereafter TW) method has been applied to a 
small, but growing, number of SB galaxies (published mea- 
surements are listed in Table 1); all cases are consistent with 
fast bars. 

Using 2-D absorption- fine spectroscopy of the SBO 
galaxy NGC 7079, Debattista & Williams (2003, in progress) 
show that the value of Qp obtained with the TW method 
is sensitive to small errors in the position angle of the disc, 
PAdisc. This raises the possibility that small errors in PAdisc 
introduce a significant scatter in TW measurements of TZ. 

Errors in PAdisc can be either simple random ones, or 
systematic ones, produced, for example, by deprojecting an 
intrinsically elliptical disc assuming it is axisymmetric. Con- 
straints on the ellipticities^ of discs come from a variety of 
measurements. The observed axes-ratios of galaxies show 



t In this paper, disc ellipticity refers to the ellipticity, ed, of the 
disc's density in its main plane. Expressions relating eo and e^, 
the ellipticity of the potential in the disc plane, can be found 
in Franx et al. (1994). Where the disc dominates the potential, 
i:d> 
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a deficit of apparently circular discs, from which one con- 
cludes that perfect oblate spheroids arc poor fits to the 
data (Binney & de Vaucouleurs 1981; Grosb0l 1985). Nev- 
ertheless, such studies find that typical ellipticities must be 
small, eo ^0.1 (MagrcUi et al. 1992; Huizinga & van Albada 
1992; Lambas et al. 1992; Fasano et al. 1993). Constraints 
on er> are improved when kinematic data are included. Rbc 
& Zaritsky (1995) defined a sample of 18 kinematically face- 
on galaxies from the Tully-Fisher relation (TuUy & Fisher 
1977, hereafter the TF relation). Using if'-band photome- 
try, they estimated typical e* = 0.051q;q2i with two arm 
spirals possibly accounting for some of this signal. Franx & 
de Zccuw (1992) showed that the small scatter in the TF 
relation requires that ej. < 0.1. Since it is highly unlikely 
that all the TF scatter is due to disc ellipticities alone, they 
concluded that a more likely limit is < eg, < 0.06. By 
analysing the residuals in the velocity-field of the gas ring 
around the SO galaxy IC 2006, Franx et al. (1994) found 
£<!. = 0.012 ± 0.026 for this galaxy. This approach has also 
been used by Schocumakcrs et al. (1997) (e$ < 0.1 for 2 
galaxies) and Beauvais & Bothun (1999), (cd ^ 0.08 for 
6 galaxies). An important uncertainty in this method is the 
viewing angle of any cUipticity. Andersen et al. (2001), there- 
fore, measured to from the discrepancies between photomet- 
ric and kinematic disc parameters of nearly face-on galaxies, 
finding an average eo = 0.05 for 7 galaxies; using the same 
method on a larger sample of 28 galaxies, Andersen & Ber- 
shady (2002) were able to fit a log-normal distribution, with 
Ini^± aine = -2.82 ± 0.73 {t5 = 0.06lo;o3). In all these 
studies, spirals may be responsible for some or all of the 
signal seen (Barnes & Sellwood 2003). Finally, in the Milky 
Way Galaxy, a variety of constraints, local and global, in- 
dependently suggest e<i) ~ 0.1, with the Sun close to the 
minor-axis of the potential (Kuijken & Tremaine 1994). 

This paper studies the effect of PAdisc errors on TW 
measurements. In Section 2 I describe the TW method and 
its main sources of uncertainty. Most of these uncertainties 
can be quantified directly from observations. However, this 
is not generally true for errors due to PAdisc uncertainties, 
so that some modelling is required. Section 3 therefore is 
devoted to setting up an A'^-body model for studying the 
impact of PAdiac errors on TW measurements. In Section 4 
I demonstrate the sensitivity of the TW method to small 
PAdisc errors and estimate the scatter in TZ expected for the 
observational level of PAdisc uncertainty. In Section 5, I con- 
sider the scatter in TZ due to non-axisyrnmctric outer discs 
on TW measurements. I also obtain a novel constraint on 
£d of early-type SB galaxies, based on the requirement that 
none of the TW measurements thus far would have found a 
value of TZ outside some range. The result is in agreement 
with previous determinations of eo for unbarred galaxies. 
In Section 6, I present my conclusions. Throughout, I pay 
particular attention to obtaining a conservative estimate of 
the scatter in 7Z due to PAdisc errors. 



2 THE TW METHOD AND ITS SOURCES OF 
ERRORS 

The TW method requires a tracer population which sat- 
isfies the continuity equation, and assumes that the time- 



Table 1. The sample of TW measurements in SB galaxies. 
The references are: Kent 1987 (K87), Merrifield & Kuijken 1995 
(MK95), Gerssen et al. 1999 (GKM99)), Debattista et al. 2002a 
(DCA02) and Aguerri et al. 2003 (ADC03). The 6 galaxies from 
Debattista ct al. (2002a) and Aguerri ct al. (2003), which have 
been analysed uniformly, constitute the ADC sample. 



Galaxy i '/'bar ^ References 



NGC 936 


41° 


66° 


1.4 ±0.3 


K87; MK95 


NGC 4596 


38° 


56° 




GKM99 


NGC 1023 


66° 


78° 




DCA02 


ESQ 139-G009 


46° 


77° 




ADC03 


IC 874 


39° 


70° 




AD COS 


NGC 1308 


36° 


60° 


o.8«:^ 


ADC03 


NGC 1440 


38° 


40° 


1 6+°-^ 


ADC03 


NGC 3412 


55° 


84° 


1 r,+0.6 
-^■'^-0.3 


ADC03 



dependence of the surface density, E, can be expressed, in 
terms of cylindrical coordinates {R, 4>) in the disc plane, as: 

S = S(J?,0-fipt). (1) 

While not all non-axisymmetric structures obviously satisfy 
the condition of equation 1 {e.g. warps), bars are well ap- 
proximated by this assumption. The TW method is then 
contained in the following expression: 

Xap = V/8mi. (2) 

Here, X = J h{Y) XT.dXdY,V = J h{Y) Mos E dX dV, i 
is the disc inclination (I use the convention i = for face-on), 
h{Y) is an arbitrary weighting function, Vios is the line-of- 
sight velocity (minus the systemic velocity) and {X, Y) arc 
galaxy-ccntcrcd coordinates measured along the disc's major 
{i.e. inclination/line-of-nodes) and minor axes, respectively. 
Equation 2 holds even when flp = flp{t), as it must, since 
the continuity equation is purely kinematic. 

Hydrodynamical studies find a narrow range in 7?- = 
1.2 ±0.2. The quoted errors and spread in TZ when measured 
with the TW method are larger (see Table 1). Important 
sources of uncertainty in TW measurements axe: 

(i) Uncertainty in Qp. To obtain Qp with the TW method, 
the most commonly used strategy is to obtain several 
absorption-lino slit spectra, for each of which V and X arc 
measured. Then plotting V versus X, one obtains Op sini as 
the slope of the best-fitting straight line. The values of X 
are usually quite well defined; however values of V tend to 
bo noisy, and arc the main source of uncertainty in Op. This 
problem can be partly alleviated by projecting slit spectra 
along the spatial direction, thereby increasing the signal-to- 
noise {S/N) ratio (Merrifield & Kuijken 1995). 

(ii) Uncertainty in Vc. Once Op is measured, Dl can be 
approximated as Vc/Op, where Vc may be assumed flat. 
However, because the tracer population must satisfy the 
continuity equation, the TW method is applied to early- 
type galaxies, which lack substantial patchy obscuring dust. 
Unfortunately, their velocity dispersions arc large, so that 
measurements of 14 require correction for the asymmetric 
drift (unless gas is present outside the bar region [Gerssen 
2002]). 
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Figure 1. The initial rotation curve of the AT-body model used. 
The dashed, dot-dashed, and dotted lines represent the bulge, 
disc and frozen halo components respectively, while the solid line 
is the full rotation curve. 








(iii) Uncertainty in as- The bar semi- major axis is some- 
times hard to measure in early-type galaxies since their baxs 
often gradually blend into the disc. The presence of massive 
bulges further complicates measurement of as. 

For concreteness, note that the mean fractional uncer- 
tainties in Q.p, Vc and an for the ADC sample (defined in 
Table 1) arc 30, 7 and 19 per cent, respectively. The resulting 
67 per cent uncertainty in TZ, averaged over all the galaxies 
of Table 1, is Atz^uuc = 0.7. (Meanwhile, the scatter of TZ for 
the full sample, which includes both an observational error 
part and an intrinsic distribution part, is A7z,obs ~ 1-0. I 
measured this value by using Monte-Carlo experiments in 
which I varied Vc and as uniformly in their error intervals, 
and varied Qp assuming its errors are Gaussian.) 

Another source of error in the TW method is errors in 
the position angle of the disc, PAdisc- Consider a slit ob- 
servation: the right-hand side of equation 2 then measures 
the flux of the tracer across the slit. However, this requires 
that the slit be exactly parallel to the X axis; for any other 
orientation, the observed velocities do not measure the full 
flux. At the same time, X, the luminosity-weighted average 
position along the slit, is rotated by the PAjisc error. The 
combination of these two effects leads to an error in the 
measured Indeed, it is surprising just how sensitive the 
TW method is to errors in PAdisc: using 2-D Fabry-Perot 
observations of NGC 7079, Dobattista & WiUiams (2003, in 
preparation) show that errors of as little as 5° in PAdisc can 
lead to errors in Qp of up to 100 per cent. Published values 
of PAdisc often have uncertainties of this order. While un- 
certainties in flp, as and Vc can be quantified directly from 
observations, errors in TZ due to PAjisc errors can only be 
modelled. 



3 MODEL AND TW MEASUREMENTS 
3.1 The AT-body system 

To quantify better the sensitivity of the TW method to er- 
rors in PAdisc, I applied it to a high resolution iV-body bar. 
In numerical simulations, Q,p can be measured accurately 



Figure 2. Logajrithmically spaced contours of the disc -|- bulge 
surface density at t = 200. The system is rotating in the counter- 
clockwise sense. 



^ tf 




Figure 3. The system after rotation through V'bar = 60°, i = 45° 
and i5pA= -|-10°. The solid line indicates the disc's true major 
axis, while the dashed lines indicate the (misaligned) slits used. 



directly from the time evolution, which makes possible a 
comparison with TW mcasuromonts at various disc and bar 
orientations and errors in PAdisc. The simulation which pro- 
duced the model of an early-type galaxy consisted of live 

disc and bulge components inside a frozen halo. The frozen 
halo was represented by a spherical logarithmic potential 



(3) 
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45°; V.^^ = 45° 



X 




Figure 4. Variation of X{Xmax) (top) and V(Xmax) (center) 
with Each line is normalized by the value of the full inte- 

gral at i5pA = 0. The bottom panel shows the resulting fra<;tional 
errors in the TW measurement of Qp using just this one slit. In 
all panels, the solid, dashed and dot-dashed lines are for i5pA = 0, 
i5pA = —5° and <5pA = +5°, respectively. The dotted vertical lines 
indicate a b ■ Other values of i/'bar and i give qualitatively similar 
results. 
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Figure 5. The TW integrals for Vbar= 45°, j = 45°. On the 

left arc shown X (top) and V (bottom) as functions of the slit 
offset. On the right, V is plotted against X, and a straight line 
fit. The solid line has slope ilpsini as measured from the time 
evolution, while the dashed line shows the best-fitting straight 
line, with slope Qxwsini. Each slit contains ^ 10^ particles; all 
errors have been enlarged by a factor of 1000 for clarity. 



where Tc is the corc-radius and vo is the asymptotic circular 
velocity. The initially axisymmetric disc was modelled by an 
exponential disc with a Gaussian thickening 



Pd{R,z) 







M 
2^' 



-R/Rd. 



R > Rt, 



(4) 



where fd is the fraction of the active mass which is in the 
disc and Rt is the radius at which the disc is truncated. 
The bulge was generated using the method of Prendergast 
& Tomer (1970), where a distribution function is integrated 
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Figure 6. As in Fig. 5 but for i5pa = —5° (top 3 panels) and <5pa 
= -|-5° (bottom 3 panels). 



iterativoly in the global potential, until convergence. For this 
application, 1 used the distribution function of a lowered, 
n = 2, polytrope, truncated at rb 

f{x,v) = C^(S) =C{[-2E{x,v)Y'^ - [-2£;^a.]i/2} .(5) 

Here C is a mass normaMzation constant and Emax = 
^tot{r'b), the total potential at rb in the disc plane. Disc 
kinematics were set up using the epicyclic approximation to 
give Toomre Q = 2.5, a value appropriate for an early- type 
disc galaxy; this leads to weak spirals, which do not interfere 
substantially with measurements of flp. Vertical equilibrium 
was obtained by integrating the vertical .Jeans equation. The 
disc and bulge were represented by 4 x lO'' equal-mass par- 
ticles, giving a mass ratio Md : Mb = /a : 1 — /a = 0.8 : 0.2. 
Further details of the setup methods used can be found in 
Debattista & Scllwood (2000). 

In units where _Rd = M = G = 1, which gives a unit of 
time {Rd/GMy^^ , the values chosen for the various param- 
eters are given in Table 2. This choice of parameters gives a 
flat rotation curve out to large radii, as shown in Fig. 1. 

The simulation was run on a 3-D cylindrical polar 
grid code (described in Sellwood & Valluri [1997]) with 
Nr xN^xNz = 60 X 64 X 225. The radial spacing of grid cells 
increases logarithmically from the center, with the outer 
edge of the grid at just over 15-Rd. The vortical spacing of 
the grid planes, Sz, was set to 0.0125J?d- 1 used Fourier terms 
up to m = 8 in the potential, which was softened with the 
standard Plummer kernel, of softening length e = 0.0125J?d. 
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t = 


Halo core radius 


rc -■ 


= 5 


Halo circulaj: velocity 




= 0.648 


Disc scaleheight 


2d 


= 0.1 


Disc truncation radius 


Rt 


= 5 


Bulge truncation radius 




= 0.78 


t = 200 


Bar semi-major axis 


a-B 


= 1.8 ±0.1 


Bar pattern speed 




= 0.296 ±0.011 


Bar speed parameter 


■R.- 


= 1.2 ±0.1 



Table 2. Parameter values of the iV-body model. 



Time integration was performed with a leapfrog integrator 
using a fixed time-step St — 0.02. 

The equilibrium set up using epicyclic theory is rather 
approximate at this high Q; nonetheless, the system quickly 
relaxes to a new equilibrium close to the initial conditions. 
The resulting axisymmetric system is unstable and forms a 
rapidly rotating bar by t = 150. Fig. 2 shows the system 
at t = 200, the time I chose for this analysis; by this time, 
the bar had gone through a period of growth and ftp had 
settled to a well defined value. The bar is strong in the disc, 
with a weaker triaxiality in the bulge. The values of the 
bar's parameters at this time are given in Table 2. Note that 
the resulting iV-body model of an SBO galaxy is reasonable, 
with a bar which is neither too weak nor too strong, having 
hb/Ra towards the upper limit of, but within, the range of 
the ADC sample. 

Since the dark matter halo is frozen, Q.p remains con- 
stant except for small oscillations produced by interference 
with weak spirals. I chose t = 200 because the spirals were 
relatively weak at this time, allowing me to measure fip with 
a minimum of interference. 



3.2 Pattern speed measurements 

For TW measurements on the A''-body system, I began with 
the disc in the xj/-plane with the bar along the x-axis, £is 
in Fig. 2. For an observer at positive z, viewing the system 
at an arbitrary orientation requires three rotations. Rotat- 
ing the system (rather than the frame), the first rotation is 
about the z-axis through an angle iphz.v, followed by a rota- 
tion about the x-axis to give an inclination i. At this point, 
the XY frame of the TW integrals is identical to the xy 
frame. A third rotation, through an angle 5pA about the z- 
axis, introduces an error in PAdisc if the observer continues 
to identify {X,Y) with {x,y). (Note that, in this definition. 
Spa > moves the assumed disc major-axis away from the 
bar's major-axis.) Prom here on, for notational convenience, 
I refer to the X and Y axes as the assumed major and minor 
axes of the system [i.e. the x and y axes), even when Spa 
/ 0. Fig. 3 shows an example of the system after such a 
series of rotations. 

I measured X and V for 0° < Vbar < 90°, 0° < i < 
90° and —90° < Spa < 90° in 11 slits covering the region 
—Ymax < Y < Ymax- Here Y,nax IS 1.2x the largest of the 
projections onto the Y-axis of the bar's 3 principal axes. This 
limited range in Y mimics the typical observational setup. 



and reduces the noise in the measurement. The values of X 
and V for each slit were obtained as: 

where 14, i and Xi are the line-of-sight velocity and X co- 
ordinate of particle i, Wi is the weight assigned to each 
particle and V — X/iesUt ™' (which corresponds to h{Y) = 
1/ J EdX, so that X and V are the luminosity- weighted av- 
erage position and velocity of each slit, as in observations). 
Except where noted, I used Wi — 1 for all particles, whether 
disc or bulge; thus V — N^ut, the number of particles in the 
slit. If A'(Xniax) and V(Xmax) represent the integrals extend- 
ing from — Xmax to Xmax, then error estimates ax and av 
were obtained by considering their maximum variation with 
XYna,x outside the bar radius. Because the number of par- 
ticles in each slit was high, these radial variations are due 
only to weak non-axisymmetric structure at large radius. In 
Fig. 4, I show A'(Xmax) and V(Xmax) for a typical slit. 

To measure the pattern speed from a set of such slits, 
I fit a straight line to V as a function of A*, as in observa- 
tions, using least-squares weights Wsut- The principal obser- 
vational uncertainty is in V and is due to photon statistics; 
I therefore used Wsnt = (cv/xZ-^sUt)"^- 

The slope of this fitted line is Qtw sin i, where I use the 
notation Qtw to distinguish from the pattern speed mea- 
sured through the time evolution. An example of such a fit 
is shown in Fig. 5, which reveals that \X\ and |V| increase 
with increasing \Y\, until they reach a maximum, and then 
decrease. Observational requirements of high S/N in mod- 
est time usually restricts slit offsets to ones at, or inside, the 
maximum in \X\ [e.g. Aguerri et al. 2003). 

I verified that the TW method accurately measures f2p 
when Spa = 0: in the range 10° < i < 80° and 10° < V'bai 
< 80°, fractional errors, jAri/fipj = KSItw — Qp)/Qp\, are 
smaller than 20 per cent, in agreement with Tremaine & 
Weinberg (1984). 

Besides this experiment, I tried various others. For ex- 
ample, in two experiments, I set Wi = and Wi = 2 for the 
bulge particles, leaving Wi = 1 for the disc ones. The results 
were consistent with those presented above, leading me to 
conclude that any plausible difference between the stellar 
mass-to-light ratio of the bulge and disc does not introduce 
large errors in Qtw- 



4 SIMPLE PAdisc ERRORS 

4.1 Sensitivity to errors in PAjisc 

Fig. 4 also plots X{Xmax) and V{Xmax) for Spa = ±5°. It 
is clear that these small errors in PAdisc change the values 
of X{Xmax) and V{Xmax) Substantially, while qualitatively 
looking similar to the Spa = case. Moreover, these changes 
are at all Xmax, particularly in the case of V{Xmax); thus, 
limiting the integrals to small Xmax does not diminish the 
error (although it does not increase it, either, unless Xmax 
is well within the bar). For this one slit, these changes gave 
an Otw which is in error by up to 100 per cent. 

In Fig. 6, I again plot the integrals as a function of 
Y, but this time for Spa = ±5°. Both \X\ and |V| reach a 
smaller (larger) maximum in the case of Spa = —5° {Spa 
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Figure 7. As in Fig. 5 but for <5pa = +5° with the disc axisym- 
metrized by shuffling the particles in azimuth. The pattern speed 
fit in the right panel is, therefore, merely an artifact. 



= +5°), while at larger offsets, the decrease in the values of 
the integrals is faster (slower) than in the 5pa = case; for 
Spa = —5°, V even switches sign. 

To begin to understand these changes, I consider an ax- 
isymmetric system. For a slit at F > 0, when 5pa = 0, the 
contribution to X and V from —X is exactly cancelled by 
that from +X. When 5pa > 0, several changes occur. First, 
+X is always closer to the galaxy center (in the disc's own 
plane), and at a smaller angle from the intrinsic major-axis, 
than is -X. Therefore \Vios{+X)\ > |Vios(-X)|, ff the rota- 
tion curve is flat, giving V a positive perturbation, which is 
further enhanced if the density profile of the disc decreases 
radially, as is generally the case. The changes in X are duo 
solely to the radial variation of the surface density; when 
this is constant everywhere, X is exactly zero at all Spa- 
Conversely, an exponential disc with small scale-length (rel- 
ative to the slit offset) gives large values of X when Spa =^ 0. 
The change in V is large already at small Xmax (sec Fig. 4), 
whereas the changes in X are more distributed over Xmax- 
This behavior is due to the fact that the integrand Y:Vios 
grows more rapidly with X than docs T,X. Indeed, for a flat 
rotation curve |Vios(-^) + Vios{—X)\ is largest at X = 0. 

Fig. 7 plots X and V at Spa = +5° for the axisymmct- 
ric disc produced by randomizing the azimuthal coordinate 
of all the particles in the iV-body model (preserving the 
average radial density profile). Even in the absence of any 
non-axisymmetric structure, misaligned slits produce non- 
zero X and V, which may plausibly be fit to a pattern speed 
where none is present. 

These extra contributions to X and V will still be 
present in the barred case, modified by the presence of the 
bar {e.g. X will still change even when the azimuthally av- 
eraged radial profile is constant, and X changes sign if the 
bar crosses the Y-axis), but fundamentally of the same char- 
acter. It is then easy to imagine that some combination of 
Fig. 5 and Fig. 7 produces the bottom panels of Fig. 6, at 
least qualitatively. For Spa = —5°, the signs of X and V in 
Fig. 7 would be reversed, which then combines with Fig. 5 
to produce something like the top panels of Fig. 6. 

Fig. 6 suggests that, when Spa < 0, it may be possible 
to recognize Spa 7^ by the large in the linear regres- 
sion. Unfortunately the most discrepant points are the ones 
at large offset; in observations, their ay will certainly be 
(fractionally) much larger than here, in which case is not 




Figure 8. Tlic maximum permitted error in PAdisc required for 
Qxw accurate to 30 per cent. Tliese have been computed in the 
range 15° < Vbar < 75°, with circles representing the averages 
and the error bars indicating the extreme cases. 



likely to be greatly increased by these points. Moreover, the 
two most discrepant points are at small \X\ and are thus un- 
likely to have been chosen for observation in the first place. 
It therefore seems likely that, in the absence of considerable 
investment in telescope time (which anyway would not catch 
Spa > 0), the error in PAdisc would go unnoticed. 

The 5° errors of Fig. 6 give errors in Qtw as large as 48 
per cent. In Fig. 8, I present the largest errors permitted to 
guarantee Qtw accurate to 30 per cent. The limits on Spa 
are quite stringent: |5pa1 < 4° is needed a,t i — 60° and the 
limit is smaller at other inclinations. (Note, however, that 
for |Af2/f2p| to be larger than 30 per cent, it is necessary, 
but not sufficient for \Spa\ to be larger than the values given 
in Fig. 8, since Spa can be either positive or negative.) 



4.2 Scatter from random PAdisc errors 

Fig. 9 plots An/fip and TZtw = {^p/^Tw)'n as functions 
of SpA. (This definition of 7?.tw ignores the errors in Vc and 
as due to Spa. These errors change 7?.tw by only a small 
amount for the inclinations of interest here.) The shaded re- 
gion in the bottom panel indicates the region of fast bars; 
it is clear that once |5pa| becomes larger than about 2°, 
values of TJ-tw scatter outside this region. Uncertainties in 
PAdisc rnust therefore also contribute to the scatter in mea- 
surements of TZ. Assuming Gaussian errors in PAdisc with 
zero mean and FWHM of 5° (2°), I found a scatter in TZtw, 
^11,6, (defined as the 67 per cent interval about the median), 
of At^^^ ~ 0.4 {^Tz^s — 0.2), as shown in Fig. 10; this is sub- 
stantially larger than the intrinsic measurement scatter at 
Spa = 0, which is only Atj ~ 0.06. Since, for the ADC 
sample, the observational root-mcan-square uncertainty in 
PAdisc is 2? 1, measurements of TZ with the TW method can- 
not directly resolve the intrinsic distribution of TZ if it is as 
narrow as hydrodynamical simulations require, even before 
other sources of scatter are considered. 

An important characteristic of the scatter is that 
TZtw < 1 may result. Since 7?. < 1 is physically impossi- 
ble (Contopoulos 1980), this may help in distinguishing the 
effects of PAdisc errors from the intrinsic distribution of TZ. 
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Ep = 0.01 = 0.05 




Figure 11. Contours of the errors in PAdisc and i resulting from assuming that an intrinsically elliptical disc is circular. The disc 
ellipticity, eo, in each case is indicated in the top-left corner of each panel. The solid contours show the errors in PAjisc, while the dotted 
contours show the errors in i. Each contour is labelled by the error it corresponds to; these are positive only for PAj;,,,, because this 
figure only considers Vdisc > 0, for the sake of simplicity. For large inclinations (near edge-on), only very small errors in PAjigc result, 
but as the disc becomes closer to face on, the errors generally become larger. The dashed lines indicate the 2 galaxies on which the 
TW method has been used with the smallest (NGC 1308) and largest (NGC 1023) apparent inclination: the bold dashed lines are the 
inclinations assumed by the corresponding authors (see Table 1), which were obtained by assuming the outer disk is circular, while the 
thin dashed lines indicate the loci of eapp = 1 — cosiapp. Where eo > 1 — cosi (e.p. NGC 1308 when eD = 0.2), the typical errors in i 
and, especially, in PAdisc become very large, up to 90°. 



5 ADDITIONAL NON-AXISYMMETRIES 

If the disc contains additional non-axisymmctric structure 
besides the bar, then this will interfere with the measure- 
ment of Q,p. If the disc non-axisymmetric density can be de- 
composed into 2 components, with different pattern speeds. 



then Otw is a luminosity and asymmetry weighted average 
of the two pattern speeds (Debattista et al. 2002b). I assume 

that the second component is a weaker non-axisymmotric 
structure and/or is at larger radius and therefore lower sur- 
face brightness, so that this type of interference will be rel- 
atively small and can be ignored. (This can also be justified 
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Figure 9. The variation of AQ/Qp (top) and TZtw (bottom) for 
small errors in PAtjiso at i = 45°. Circles, triangles and filled 
squares are for Vbar = 30°, 45° and 60° respectively. The dotted 
lines in the top panel represent errors of 20 per cent, while in 
the bottom panel, the shaded region indicates 1.0 < TZtw < 
The crossed open circles in the top panel represent a system with 
i and Vbar as in NGC 7079. 
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Figure 10. The distribution of TZtw for various distributions 
of random errors in PAdigc- The solid line shows the distribution 
without PAtjisc errors, while the dot-dashed and dashed lines show 
the distributions resulting from Gaussian errors of zero mean and 
FWHM = 2° and 5°, respectively. The dotted Une is also for 
Gaussian errors with FWHM = 5°, but uses = <tv~^ to 

measure Htw- Ea<;h line has been rescaled vertically for clarity. 
The distributions represent averages over 30° < i < 70° and 
10° < i/ibar < 80° , and are not substantially changed by modest 
changes to these limits. 



by noting that the weak spiral structure at large radius in 
the A'^-body model does not introduce substantial errors in 
ftxw.) Instead, I concentrate only on the effect these sec- 
ondary non-axisymmetric structures have on Qtw due to 
the errors they introduce in the measurement of PAjisc. 



5.1 Elliptical discs 

In all cases in which the TW method has been used, PAdisc 
has been measured from surface photometry under the as- 
sumption that the disc is intrinsically circular. When the 
disc is elliptical, deprojectirig with this assumption gives rise 
to errors in i and PAdisc, as shown in Fig. 11. These errors 
lead to further scatter in TZtw- 

To study this scatter, I assumed that, at large radii, en 
and Vdisc (where i/idisc is the angle of the elliptical disc in 
the plane of the disc relative to the line-of-nodcs) arc both 
constant, and computed the apparent PAdisc (PAopp) and 
apparent i (iapp) resulting from the assumption of a circular 
disc. I used these to measure the apparent circular velocity 
{Vc,app) and bar semi-major axis {aB,app)- I then obtained 



f^Tw sin iapp as the slope of the best-fitting line to {X,V), 
from which I measured 72.tw = Vc,app/{aB,app ^tw). By 
assuming that the bar is infinitely narrow, I measured the 
apparent bar PA in the disc plane, 4'b,app, and then averaged 
7^Tw over 30° < iapp < 70°, 10° < i^t.app < 80° and -90° < 
V'disc < 90° . Fig. 12 plots the resulting distributions of 7?.tw 
obtained for various constant eo- The ellipticity-induced 
scatter, A-jz^^, grows rapidly with eo [Atz^^ ~ 0.2, 0.6 and 
0.9 for eo = 0.01, 0.05 and 0.1 respectively), with most mea- 
surements of 72.TW outside the range 1.0 < TZtw < 1-4 once 
eD = 0.1. The distinctive peak to 7?.tw < 1 for the larger 
values of eo is due to the fact that the distribution of Spa, at 
fixed iapp, has peaks near max(|(5pA|). The peak at TZtw < 1 
is higher than that at 7?-tw > 1 because TZ cx fJp ^. 

Fig. 12 also shows the distribution of 7?.tw resulting 
from the ellipticity distribution of Andersen & Bershady 
(2002) for later-type unbarred galaxies. The two largest val- 
ues of €d in their sample of 28 were en = 0.232to'o64 and eu 
= 0.165 ± 0.083 (Andersen 2002, private communication). 
As can be seen in Fig. 11, large values of eo produce PAdisc 
errors as large as 90° in the (apparent) inclination range 
of interest, which would result in very large errors in Qtw- 
Therefore I truncated their distribution at eo = 0.1 and 0.15. 
The resulting scatter is A-r,,^ ~ 0.5 and ~ 0.6 respectively. 

To compute an upper limit for the characteristic eo of 
SBO galaxies, I define Pf as the probability that all mea- 
surements will result in 0.5 < 7?-tw < 2.5, a range outside 
which, at the 67 per cent interval, none of the measure- 
ments of Table 1 fall. Then, for that sample, I compute Pf 
by matching iapp and ipb,app to the observed values and av- 
eraging over V'disc, obtaining Fig. 13. The probability of hav- 
ing found 7?-TW loss than 0.5 or greater than 2.5 for one or 
more of these galaxies exceeds 90 per cent (75 per cent for 
TZtw > 5.0) if eo > 0.07 for all of them. (The strongest 
constraints come from the low inclination galaxies, while 
NGC 1023, which has the largest inclination of this sample, 
does not constrain eo at all, up to 0.1.) This upper limit 
on the disc ellipticity is in rough agreement with previous 
measurements {e.g. Franx & de Zeeuw 1992) for unbarred 
galaxies. 
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Figure 12. The distribution of 7?.tw resulting from elliptical 
discs. The thin solid line shows the intrinsic distribution when 

eo = 0, while the thin dashed, dot-dashed and dotted lines show 
the distributions resulting from errors caused by en = 0.01, 0.05 
and 0.1, respectively. The thick solid and dashed lines are based 
on the empirical distribution of Andersen & Bershady, with a 
maximum cd of 0.1 and 0.15 respectively. Each line has been 
rescaled vertically for clarity. The distributions represent averages 
over 30° < iapp < 70° and 10° < ipb,app < 80°, and are not 
substantially changed by modest changes to these limits. 
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5.2 Rings 

In Section 5.1, I assumed that ^disc is uncorrclated with 
tphai- Correlations between i/'disc and i/'bar niay be introduced 
by the outer rings often seen in SB galaxies. Two main types 
of outer rings are possible {e.g. Buta 1995): Ri, which are 
aligned perpendicular to the bar, and R2, which line up with 
the bar. Galaxies selected for TW measurement do not con- 
tain strong rings, but conceivably weak rings might have 
been overlooked. To consider their effect on TW measure- 
ments, I simply set -0bar = '/'disc (for rings of type R2) and 
V'bar ~ ^disc+90° (for rlngs of type Ri) and proceeded as 
for Fig. 12. The results, unsurprisingly, sfiowed that rings of 
type R2, which lead to 5pa < 0, produce 7?.tw < T^, while 
rings of type Ri lead to TZtw > TZ- Buta (1995) found mean 
CD of 0.26 and 0.13 for rings of type Ri and R2 respectively. 
If such rings had been present in the sample of Table 1 , then 
the scatter in 7?.tw would have been significantly higher. 



5.3 Spirals 

Recently, Barnes & Sellwood (2003) have questioned the in- 
terpretation of discrepancies between photometric and kine- 
matic inclinations and PA's as resulting from disc elliptic- 
ities. Instead, they found evidence that spirals, or similar 
non-axisymmetries, produce these discrepancies. They re- 
ported an average PAdisc uncertainty of about 4° for earlier- 
type galaxies. 

PAdisc errors of this type will produce scatter in 7?.tw 
in much the same way as do random PAjisc errors. For a ~ 
4°, I found a resulting scatter Ajz^spr — 0.7. However, the 
sample of galaxies used by Barnes & Sellwood (taken from 
Palunas & Williams [2000]), excluded galaxies as early as 
SO, so this value is somewhat uncertain and is probably an 
over-estimate. 



5.4 Warps 

While most disc galaxies arc coplanar inside R25 (Briggs 
1990), examples of warps inside this radius are not unknown. 
One extreme case is the interacting galaxy NGC 3718, which 
has a warp of about 80° at R25 (Schwarz 1985). However, 
such strongly interacting galaxies are usually not selected 
for TW studies. Furthermore, the large velocity dispersions 
of early-type gala:xies serve to stiffen their stellar discs (De- 
battista & Sellwood 1999), so that any warps inside R25 are 
generally small. Therefore warps probably do not introduce 
significant scatter in TW measurements. 



Figure 13. The bottom panel plots the probability, Pf, that none 
of the 8 galaxies of Table 1 is outside the range 0.5 < T^tw < 2.5 
for £□ fixed for all galaxies. The top 8 panels show the distri- 
butions of T^TW (filled histograms and bottom scale) and 5pa 
(open histograms, top scale) produced by matching ipb^app 
iapp for the galaxies in Table 1 under the assumption that ed 
= 0.07 (where Pf < 0.1). 



6 DISCUSSION AND CONCLUSIONS 

6.1 How realistic are the error estimates? 

How realistic are these estimates of AQ/^lp and the At^'s? 
Since, for Spa < 0, the values of X and V are not all close to a 
straight line (see Fig. 6), a poor choice of M^siit could lead to 
excess scatter. At 6pa = 0, I obtained the smallest |An/f2p| 
with VKsiit ~ o-v~^ , which is defined only from the variations 
of V with Xmax- This is unsurprising, since av represents the 
full uncertainty in V. All other definitions of Wsnt produced 
larger errors. In particular, while Wsxit = crv~^ gives a mean 
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AQ/^lp of 3 per cent, Went = (o"v/\/-^siit)~^ produces a 
mean AQ/Qp of 7 per cent. 

However, when Spa 7^ 0, VFsiit favoring slits with small 
offset, which generally acquire fractionally smaller pertur- 
bations, produces smaller scatter. Fig. 10 compares the dis- 
tributions of 7?.Tw from random Gaussian PAdisc errors of 
FWHM = 5° as obtained using W^t = {crv / VKiit)~^ 
and Wsut = crv~^- The former produces a smaller scatter, 
due mostly to the reduced noise at TZtw > 71, i-e. at Spa 
< 0. I tried other definitions of VI4iit, including ax^^, equal 
weights, V, and various combinations of these. I also tried 
using only 3 slits (the central one and either the two with 
the largest lA"] or the two flanking slits), as is often done in 
observations. These always gave larger scatter, typically by 
20 per cent or more. I therefore used Wsut ~ {o'v / VNsm)^'^ 
everywhere in this paper to compute Af2/f2p and the Ar's. 
Thus I am assured of a conservative estimate of the scatter, 
while also matching better the main source of noise in the 
observations: the photon statistics. 

Since I have used only one simulation to estimate the 
scatter, I need to show that this simulation does not over- 
estimate the errors in VItw that real galaxies would suffer. 
Perhaps the most important parameter affecting the size of 
the scatter in Htw is as/Rd, as described in Section 2. A 
series of experiments with razor-thin, flat rotation curve, ax- 
isymmetric exponential discs showed that, indeed, the scat- 
ter in f^Tw due to random PAaisc errors increases as -Rd 
decreases. Since my model SBO has a value of as/Rd that 
is towards the upper end of those in the ADC sample, my 
measurements of A-n^s and A-jz^e probably underestimate 
somewhat the scatter which the same PAdisc errors would 
produce in real galaxies. The same conclusion resulted from 
a test with a lower quality (102K particle) simulation hav- 
ing a larger bar {as/Rd = 2.6); for random Gaussian errors 
of FWHM = 5°, this bar produced A-R,i = 0.3 versus 0.4 
for the shorter bar used in this paper. 

The trend with Spa seen in Fig. 6 is in the same 
sense as was found by Debattista & Williams (2002) for 
NGC 7079. Fig. 9 plots Afi/fip for the same projection as 
NGC 7079. The errors in Qtw due to Spa for NGC 7079 
{as/Rd = 1.5 ±0.2) reported by Debattista & Williams are 
perhaps a little larger than those computed here. Gratify- 
ingly, the error estimates produced by the iV-body model 
are not unrealisticly large. 

6.2 The ellipticity of early-type barred galaxies 

The ellipticities of SO galaxies are poorly constrained. From 

photometry only, Fasano et al. (1993) found that they could 
not rule out that they are perfectly oblate. The two SO galax- 
ies with directly measured ellipticities, IC 2006 (Franx et al. 
1994) and NGC 7742 (Rix & Zaritsky 1995) both have small, 
possibly zero, ellipticity (e* = 0.012 ± 0.026 and 0.02 ± 0.01 
respectively). The ellipticities of SB galaxies arc not irmch 
better constrained, undoubtedly because they require a dis- 
tinction between the inner, bar-dominated, region and the 
outer parts. Photometry alone, therefore, is of limited use, 
and kinematics also are needed. Unfortunately, most TF 
studies have avoided SB galaxies. Debattista & Sellwood 
(2000) showed that the small fraction of bright {Mi < -21) 
SB galaxies contaminating the sample of Mathewson & Ford 
(1996), who selected against SB galaxies, satisfies the same 



TF relation, and has the same scatter, as the unbarred 
(SA) galaxies. Sakai et al. (2000) cahbrated the TF rela- 
tion of nearby galaxies with Cepheid distances; their sample 
of 21 galaxies contained a more representative fraction of SB 
galaxies, at ~ 30 per cent. The resulting TF relation, includ- 
ing the scatter, also was identical for SA and SB gala^xies. 
Thus we may suppose that the TF-based constraint of Franx 
& de Zeeuw (1992), e# < 0.1, also holds for SB galaxies. 

The constraint obtained here, eo ^ 0.07, is in rough 
agreement with the constraints for SA galaxies. However, 
an important possible bias needs to be pointed out. The 
ADC sample of 6 galaxies explicitly excluded galaxies for 
which, at large radius, the observed PAdi^c changes substan- 
tially with radius. From a sample of 11 galaxies for which 
they obtained surface photometry, one (Aguerri 2002, pri- 
vate communication) was excluded for this reason. If either 
eD or V'disc changes with radius, then the observed changes 
in PAdisc will typically be greater in galaxies with larger 
mean to- Thus the cut on the size of PAdisc variations may 
have introduced a bias in the ellipticity distribution of the 
ADC sample; on the other hand, large variations in PAdisc 
may have been caused instead by spirals or by a warp. 

Although these constraints on SB galaxy ellipticities are 
consistent with the constraints on SA galaxy ellipticities, 
this does not mean that their ellipticity distributions arc 
the same, since both the TF and the TW constraint obtain 
only upper limits on eo- 

6.3 The intrinsic distribution of TZ 

Hydrodynamical simulations of SB galaxies find a narrow 
range of TZ — 1.2 ± 0.2. The presently observed distribution 
of 7?.TW is dominated by the observational uncertainties in 
f^Tw, os and Vc- Nevertheless, it is clear that all 8 galax- 
ies measured so far are consistent with the range found in 
hydrodynamical simulations. In their A'^-body simulations 
with cosmologically motivated initial conditions, Valcnzuela 
& Klypin (2002) found that bars with 7?. = 1.7 were pro- 
duced, which they considered to be consistent with the ob- 
servations. Indeed, for 4 of the 8 galaxies listed in Table 1, 
TZ — 1.7 is within the error interval. However, 3 of these 4 
gala^xies are the ones with the largest error bars, and the 
fourth galaxy is only just barely consistent with this value. 
For the ADC sample, which have well-determined PAdisc 
uncertainties, the root-mean-square uncertainty in PAdisc is 
2?1. From the results of Section 2, the corresponding scatter 
in TZtw, excluding any contribution due to disc ellipticity, 
should be Afi,s — 0.4. Allowing for this scatter, it seems 
possible that 1.7 is outside the intrinsic range of TZ. 

For a crude estimate of the intrinsic range of TZ, sup- 
pose we can write A^^^^^ = A|j_i„t + A^^^g + A^^^ + A^^^„^, 
where A-jz^obs is the observed scatter, A-jz^int is the intrinsic 
range of TZ, A-n^g is the scatter due to random PAdisc er- 
rors, Atj.c is the scatter due to disc ellipticity and A-jz. unc 
is the scatter induced by uncertainties in the measurements 
of flp, as and Vc. All these A-r,'s are assumed to be 67 per 
cent intervals. (Other sources of scatter, such as direct in- 
terference from spiral or other structure, small errors in slit 
orientation, etc., may be present but arc assumed here to 
be unimportant.) From Section 2 I get that A-iz^obs — 1-0 
and ATz,unc — 0.7, while from Section 4 I get A-jz^s — 0.4. If 
£d = for all galaxies, then ATz,int — 0.6, while the distri- 



© 0000 RAS, MNRAS 000, 000-000 



On Position Angle Errors in the Tremaine- Weinberg Method 11 



bution of Ed of Andersen & Bcrshady (2002), truncated at 
£d = 0.1, produces A7j.,i„t — 0.3. If, on the other hand, the 
interpretation of Barnes & Sellwood (2003) is correct, then 
A7j,e = 0, but it is replaced by A ^ 0.7. It therefore 
seems possible that the intrinsic range of TZ for early-type 
gala^Kies spans a range similar to the later-type gala^xies. 

Unfortunately, the sample size is still too small for a 
proper statistical test of this suggestion. If correct, then the 
fact that SB galaxies have the same distribution of TZ as the 
more gas-rich later-type SB galaxies requires that gas is not 
dynamically very important for the evolution of Q.p. 



6.4 Future work and conclusions 

The current sample of TW measurements is still quite small, 
so it is not unlikely that, in the future, more nicasuromonts 
will be obtained. The results of this paper can be read as an 
endorsement of careful surface photometry of target galax- 
ies to accurately measure PAaisc. Inclinations in the range 
50° < i < 60° are preferable, since they are less sensitive 
to errors in PAdisc- For statistical studies, especially to con- 
strain the distribution of TZ, it would be very useful if future 
studies were to report their uncertainty in PAdisc. Galaxies 
with strong outer rings do not make good candidates for TW 
nicasurcrncnt because of the inherent uncertainty in PAdisc, 
and should be avoided. If the TW method is ever to be used 
on late- type galaxies, perhaps in the infra-red {e.g. Baker 
et al. 2001), care must be taken that the presence of spirals 
does not lead to excess errors in PAdisc. 

The pattern speed of triaxial elliptical galaxies is a mat- 
ter of theoretical speculation. Because of the large velocity 
dispersions and low stellar streaming velocities, it is gener- 
ally thought that their pattern speeds must be small. Mea- 
surement of their pattern speeds would be very interesting, 
but unfortunately, application of the TW method to ellipti- 
cal galaxies is likely to be accompanied by significant uncer- 
tainty in their intrinsic orientations (amongst other difficul- 
ties). Thus TW measurements of their pattern speeds may 
have large uncertainties. 

I have shown that errors in PAdisc lead to significant 
error in TW measurements. For the observational level 
of random Gaussian errors, the resulting scatter in TZ is 
1^71,6 — 0.4. If barred galaxies are also modestly elliptical, 
then the total scatter increases further, depending on the 
distribution of cd • Given the observed range of TZ, this sug- 
gests, therefore, that the gas-poor early-type galaxies have a 
narrow distribution oiTZ ^ 1.0 — 1.4, not much different from 
gas-rich late-type galaxies, as determined by independent 
means. This result would imply that gas is not dynamically 
important for the evolution of bar pattern speeds. 
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